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Abstract 

New version of the MSSM scales is discussed. In this version \x <C Msusy ~ 
Mq ~ M1/2, where /i is the Higgsino mass, Mq is the mass scale of sleptons and 
squarks, M 1 / 2 is the mass scale of gaugino. Renormalization group motivation of 
this MSSM version is proposed. The radiation corrections give main contribution to 
the splitting of neutralino and chargino masses if Msusy > W 7 GeV. Calculation 
of the mass difference has led to the value M ± — M^o ~ 360 MeV at M v ~ u ~ 

Xi Xi,2 X r~ 

3TeV. The mean life time of charged Higgsino in this model is r ± = 0.4 x 10~ 9 s. 

Xi 

The analysis of the residual neutralino concentration in the cosmological plasma is 
developed and the value of Higgsino mass is estimated: M x ~ 3TeV. Formation 
of the residual neutralino concentration occurs in the high symmetric phase of the 
cosmological plasma. Problems of the relic Higgsinos searches at the underground 
laboratory NUSEL and Galactic neutralino annihilation radiation at the satellite 
detector GLAST are discussed. It is shown that the neutralino-nucleon cross section is 
spin-dependent. The annihilation of Galactic neutralinos occurs in Z°Z° and W + W~ 
pairs in t- and s-channels, and also in lepton-antilepton and quark- antiquark pairs in 
the s-channel. The contribution to the annihilation spectrum from first two processes 
is calculated with using experimental data about hadron multiplicities at the mass 
surfaces of Z°, VK^-bosons; the calculation of quark-antiquark s-channel contribution 
is based on the phenomenological model of hadron multiplicities at ^/s = 2M X with 
using negative binomial distribution. 



1 Introduction. 

As it is known, a theory doesn't fix the MSSM scale hierarchy. Both neutralino and chargino 
have different properties depending on this hierarchy, therefore predictions for experiments 
are also different. Neutralino is nearly bino in generally accepted variant of the theory, and 
neutralino interactions with matter mainly occur through the scalar quarks' exchange since 
their masses are close to the neutralino mass. The splitting between neutralino and chargino 
masses is large enough and its value has order a few GeV. It allows to register these particles 
through the lepton and hadron cascade decays. The standard MSSM version implies an 
existence of the Supergauge Desert between electroweak and GUT scales. But another 
variants are also possible. In particular, we are interested in in-depth study of the MSSM 
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variant, where the mass scales of squarks and sleptons and the gaugino scale were arranged 
far from electroweak scale in multi-TeV region (so-called the Split Supersymmetry) pQ, [2|. 
Formally Higgsino mass remains arbitrary so it may be positioned at the EW scale. It's so- 
called the light Higgsino scenario. We study in details the properties of this model, which is 
interesting because in the case ideas of SUSY are conserved. It should be especially actual 
if LHC results for the SUSY effects' observation will be negative. 

2 Objects of investigation and its properties. 

The experimental testing of the supersymmetry hypothesis is one of the general purposes 
of Large Hadron Collider [3], 0j, [Sj. Main object of investigation is the Minimal Su- 
persymmetric Standard Model, which represents the supersymmetric generalization of the 
Standard Model. Experimental program of supersymmetry searches on the LHC is based 
on the SUSY breaking scale estimation Msusy ~ 1 — 2TeV. There are four theoretical 
arguments that have lead to the estimation: 1) the MSSM is the simplest supersymmet- 
ric model compatible with well-known experimental data; 2) the mass spectrum of Higgs 
bosons is stable; 3) precise convergence of invariant charges at one point; 4) in the frame- 
work of MSSM we have natural interpretation of Dark Matter. 

As it is known, after spontaneous breaking of the gauge symmetry supersymmetric 
partners of Higgs and electroweak gauge fields forms two Dirac electrically charged particles 

- charginos X121 an d four Majorana electrically neutral particles - neutralinos Xa> a = 
1, 2, 3, 4. Lightest supersymmetric particle Xi is main candidate to be the cold Dark Matter 
component [H], 0. Generally Xi is the superposition of U(l) gaugino B ("bino"), 
neutral SU(2) gaugino W3 ("wino") and two Higgsinos h®, h®. Structure and properties of 
chargino and neutralino depend on relations between characteristic MSSM scales. Strong 
theoretical arguments in favor of one or another hierarchy of scales are absent. 

In the next section will be shown that only two variant of MSSM are theoretically 
natural. First variant - the light gaugino scenario (bino or wino or their mixture), in which 

\fi\ > M ~ Mi/a ~ Msusy > M EW , (1) 

- is well investigated early [TU], |Hj, [12], [IS]- 

We consider second alternative variant of MSSM, in which the Higgsino mass scale is 
the nearest for EW scale, but Msusy is in the multi-TeV region - so-called light (separated) 
Higgsino scenario: 

M ~ Mi/ 2 ~ Msusy > M > M EW . (2) 

3 Renormgroup analysis. 

In our model we shown precise convergence of invariant charges at the scale compatible with 
proton mean life. In the known MSSM scheme this convergence of invariant charges at the 
Mqut provides the MSSM building-in to the Supergravity and super string theory. In the 
case the scale of SUSY breaking is no so large: M S usy ~ Mi/ 2 , so quadratic divergencies 
cancel near M#. 

We take into account all degrees of freedom in the SUsusy(5) except of quark-lepton 
bosons. It have been found that the states, which are near Mqut, are very important for 
the renormgroup analysis too. 
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a^(M z ) = -a-\M z )(l - sin 2 9 W (M Z )) = 59.0132 =f (0.0384) sin 2 ew . ± (0.0124) a , 



Initial values of running constants are fixed at the M^-scale: 
oT\M z ) = 127.922±0.027, a s (M z ) = 0.1200±0.0028, sin 2 9 W [M Z ) = 0.23113±0.00015. 

To solve renormgroup equations the following initial values are used: 

3 
— ( 

5 

a 2 \M z ) = a~ 1 {M z ) sin 2 9 W (M Z ) = 29.5666 ± (0.0192) sin2(?w . ± (0.0062) Q (3) 

a" X (M Z ) = a-\M z ) = 8.3333 ± 0.1944. 

In the one-loop approximation invariant charges at the scale Q 2 are linked with its' 
values at the scale Qi. 

oci\Q2) = *r\ Ql ) + |ll n ^, bl = 5>,, (4) 

3 

Here all states with masses Mj < Q2/2 at Q2 > Qi are summed. 

Extra degrees of freedom are taken into account in our analysis, namely: singlet quarks 
and their superpartners (D L , L>l), (Dr, D r ) from superhiggs quintets of SUsusy(5); gauge 
superfield $l) from 24-plet in adjoint representation of SU{2) and chiral superfield 
in the adjoint representation of SU(3). In the minimal SUsusy{5) we have 
M 5 = (Md, M d ), M 24 = (M$, M^,, M$, M|,) and all of them are generated by interaction 
with Higgs condensate at Mgut- Inequality M 5 , M 24 < Mgut, in principle, can be provided 
with precision in 1 — 2 orders. 

So, one-loop invariant charges at the scale q 2 = {2M G ut) 2 depend on all characteristic 
MSSM scales: 
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Here Mo = (M3, Mr) is the "common" mass of scalar quarks and leptons, averaged over 
chiralities and generations; M t is the t-quark mass; others mass parameters were introduced 
above. In compliance with internal feature of the MSSM in (jSJ) it supposes that the lightest 
Higgs boson mass is close to M z , but masses of others Higgs bosons H, A, H ± lie at the 
M H scale. Expressions (j5J) do not depend on placing of M 24 , M 5 , M , M1/2, /i, M H at the 
energy scale. If masses of singlet superquarks and residual Higgs superfields are equaled to 
Mgut they are eliminated from renormgroup equations. 
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At every scale corresponding degrees of freedom are linked and the equaling of all charges 
at Mqu T leads to equations: 

\ 2/9 / \ 1/3 

M' 1/2 ) " z \M' 1/2 



where 



U<» i = Ak { \! [ — | . y.= Bk 2 M z ] . (G) 

1/12 , v 1/3 



L _ 7^-1/12^1/3 - ( M l\ ( M GUT\ 

1 — /v sf ^GUTl 



\ m qJ \ 1V1 GUT J 

k^k^ki'ik-^ - ( ^v /4 r^v /4 r^r & 
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,t ^ gut2 y Mg j y M j y M ^j y McuT 
M' 1/2 = (M^Mg) 1 / 2 , M' GUT = {M^M^iM^M^ < M GUT , 

M GUT = (M | M$) 4/3 < M OUT, 

( 8 ) 

A = exp (^-^ a -\M z ) - 3a 2 \M z ) - 2a^(M z )) - ^ln2j = (1.57xJ;|») ■ 10 14 , 

B = exp l-(5a^(M z ) - l5a 2 \M z ) + 7a^(M z )) + — ln2j = (2.0x°;*jj) ■ 10 3 . 

The mass ratios that are labelled as K are quantities, which are larger than unity. 
All parameters Kquti, K GUT2 both in the MSSM, and in the SUsusy{5) are not under 
the theoretical control. Due to experimental restriction M H > 315 GeV it is supposed 
that 1 < Kquti — Kqut2 < 10 and Km lies in the region 2 < Km < 10. Parameters 
K~i, K~yy are fixed by renormgroup evolution from Mqut to M , Mx/ 2 . We suppose that 
l.5<K~^K~ gW <2.5. 

The studying of M[, 2 , \i as functions of M GUT is the aim of analysis: 

M[ /2 {M GUT ) = {Ak^Mf 2 x -1 ^{Mqut) = — ^ /2 x M^ 2 T . (9) 

M GUT \Aki) I M z 

We used known restrictions for the proton mean life and for Msusy'- ( r p > 10 32 years 

/2 



at Mqut > 10 15 GeV); (M 5[/5y ~ M' > 100 GeV at M GUT < 3 • 10 16 GW) 



i/ 

From the analysis two variants of parameters (jHJ) were found. In the first case \x <C M[/ 2 , 
in the second one // ^> M[\ 2 that is shown at Fig. [T] 



4 Analysis of mass spectrum. 

Formally the Mjj scale may be arbitrary. In this scenario light particles Xi 2 anc ^ xf have 
a structure of Higgsino with masses 



M|(l + sign(n) sin 2/?) , ^ cos # sin 6 1 \ M 
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M ± ~ - ^£ ( J^L + s^n^) sin 2/3 J « |/i|. 




Figure 1: Two variants of the MSSM scales - bino-like LSP \i S> M[, 2 and Higgsino-like LSP 
M«M{ /2 . 

Heavy particles Xs, 4 and xf have a structure of gaugino, their masses lies at M SUSY scale: 

M x o^M § , M x o^M Wl M xt ^M w . (11) 

Interactions through virtual gaugino, squarks and sleptons and heavy Higgs are sup- 
pressed by large masses in denominators. Basic physical effects describing by the La- 
grangian of light Higgsino Xi 2 > Xi interactions with photons and vector bosons: 

AL = (eA, - - 2sin 2 0)<) x^YXi + i^-^M^xl + ^7 M X?)+ 

V L COS o / Z COS v 

(12) 

+ xS)Vxr + ^w-xmxl + xg). 

Tree constraint on the mass splitting between chargino and neutralino M ± — M o > m e 
with (jlUJl gives some constraint on supersymmetry breaking scale: Msusy < 10 7 GeV. But 
if Msusy > 10 7 GeV the radiative corrections dominate in formation of the neutralino- 
chargino mass spectrum. Combination of all mass operators in the one-loop approximation 
may lead to the finite integral on the mass surface: 

ie 2 M 2 7 f (q- M y )dq 

AM V = - z 1 w x> H 



87T 

For M x > M z we have 



q%q* - Mf)[(g + pf - Af|] 
a(M x )M z 



AM X ~ v ^ (13) 



To obtain the value a at M x scale by renormgroup methods we use estimation Mh 
M t ~ 200 GeV UU, [T5| and M x ~ 3TeV that will be obtained below. We found AM 
360 MeV. The calculation of chargino mean life time gives 



x — 



V = ^(AM x )- 5 = 0.4xlO- 9 ,. (14) 

If we refuse from Supergauge Desert we haven't theoretical considerations on the sep- 
arated Higgsino mass. Nevertheless we can attempt to make an important step to solve 
this problem when we use the hypothesis that neutralino is a carrier of Dark Matter in the 
Universe. 



5 Cosmological estimation of the separated neutralino 
mass. 

In discussing versus of the theory physics of neutral and charged Higgsinos fixes by the 
Lagrangian (fT2|) and the mass spectrum (|10jl with the radiation splitting (|13jl . Therefore 
any quantitative predictions of the theory depend on two parameters /i, M# only. Concrete 
value Mh isn't essential. The basic parameter of the model with separated Higgsino is [l. 
Let imply that neutralino is a carrier of Dark Matter in the Universe. Irreversible neutralino 
annihilation starts in the cosmological plasma at the moment to and the temperature To, 
when mean energy of relativistic quarks and leptons compares with neutralino mass: if ~ 
3T = M x . Residual density of neutralino at t t describes by the asymptotic solution 
of evolutional equations: [H] , [7J , [TB] : 



-dT + 3Hn *--2"> 



x + 3Hn x = --n x (av) ann , 

(15) 



3H* = 8nGwT\ H = i ■ ^, wT* = 

a dt a 4 

Here (av ) ann - kinetic cross section of the annihilation, which approximate to the constant 
in our model; G - gravitational constant; w = w(T) - statistical weight of plasma at the 
annihilation epoch. With all degrees of freedom in the two-doublet Standard Model we have 
w = 4437r 2 /120. This value was used us for analysis of annihilation at the temperature 
exceeding the characteristic scale of Standard Model. 

The asymptotic solution of equations ()15j) at times, exceeding the moment when anni- 
hilation terminates, is 

xl ' V <>(«) / V 3 / (H.™ 

where T ~ M x /3 is the temperature at the moment of annihilation beginning. Annihilation 
epoch begins at the moment 

1 / 3 \ 1/2 

(17) 

and ends when T\ ~ M x /20. To obtain the value of mass density of the neutralino gas in 
the present-day Universe p x = M x n x we use the standard data for cosmological neutrino. 
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Evolution of the neutrino component in the cosmological plasma may be considered as adia- 
batic since the energy isolation occurs mainly in non-neutrino channels of the hadronization 
and annihilation processes at t > to- Therefore the ratio of scale factors in (|T6|) we can 
replace by ratio of neutrino gas temperatures: 

^4 = % T - (t?I T - = L676 x 10 ~ 13 G ^ 

a[t v ) T V 11 / 

Here t v is the age of the Universe; T u = T„(t^), T 7 is the temperature of the relic gamma 
radiation. 

Finally for density of stable neutralinos in the epoch, characterizing by the temperature 
of relic neutrino T„, we obtain: 

M^(^Y /2 ^, (18) 

lo \ o J {<JV) ann 

From the recent WMAP data J7j, JH] the mass density of Dark Matter in the present- 
day Universe is equal: 

p DM (t v ) ~ (0.23 ± 0.04)p c ~ (0.94 ± 0.34) x 10~ 47 GeV 4 . (19) 

Here p c = (4.1 ± 0.8) x 10~ 47 GeV 4 is the critical density of the Universe. An assumption 
about the Dark Matter in the Universe is the neutralino gas provides in 

Pxitu) = p DM {tu). (20) 

To estimate the value of neutralino mass M x from (J20j) we should obtain the expression 
for kinetic cross section of neutralino annihilation in the cosmological plasma [av ) ann as 
a function of M x . The problem is that we don't know in which phase of the cosmological 
plasma neutralino annihilation occurs. We analyze two scenario consistently. 

1. Annihilation in the low symmetric phase. Formation of the residual neutralino 
concentration occurs in the low symmetric phase (LS-phase) of the cosmological plasma 
that is after the electroweak transition if neutralino mass M x < 3T e vy ~ 300 GeV . In 
standard scenarios of the cosmological evolution of neutralino as gaugino in analysis of 
the Majorana neutralino annihilation to massless fermions for execution of the Pauli's prin- 
ciple we must take into account of thermal corrections as a result the full annihilation cross 
section is a function of the temperature. Formally we have two variants: the annihilation 
of short- living particles (r x < t ) and the coannihilation of long-living particles (r x ^> to) 
[T2] . The carried out analysis shown that for existing of the separated Higgsino scenario 
we must forbid the s-channel annihilation into quarks and leptons that is to make second 
Majorana neutralino X2 unstable with mean life time is smaller than 10~ 10 s. It is possible 
if Msusy scale doesn't exceed 10 5 GeV. But this scenario assumes the fine tuning of the 
neutralino mass to W-boson mass. Therefore 

M- ~ M Y o ~ M w , 

that is contrary to well known experimental constraints |19j . Formally equations for an- 
nihilation in the LS-phase contain the second solution corresponding to neutralino mass 
in order a few TeV. But such scale is incompatible with the assumption about neutralino 
annihilation in the LS-phase of the cosmological plasma. So this assumption doesn't reduce 
to self-consistent result. 
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2. Annihilation in the high symmetric phase. As it is known, at temperatures, 
exceeding the temperature of electroweak transition T > T e w ~ 100 GeV, plasma is in the 
high symmetric phase (HS-phase), in which Higgs condensate is absent. The characteristic 
feature of the HS-phase is the massless of all particles (more precisely, their masses is 
m < T) except Higgsino. Physical states here is quants of gauge fields B, W a (a = 1, 2, 3) 
and chiral fermions. The mass splitting disappear in the Higgsino family. As a result 
neutralino and chargino degrees of freedom incorporate to the fundamental representation 
of the 577(2) group that is the unified Dirac field x> m which internal states classify over 
quantum numbers of the restored SU(2) symmetry and they are dynamic equivalent. At 
the HS-phase instead of (fT2^) we should use the Lagrangian 

AL X = ^B^rx + \g2W;xYr aX , (21) 

which adds to the Standard Model Lagrangian written in terms of gauge and chiral fields: 

1- _ 1 2 1 - 

L S m = --giB^l L ^l L - giB^e R ^e R + -giB^q L ^q L + -giB^u R ^u R - -g\B^d R ^d R 

+ \g2Wft L Yr a l L + \g2W«q Ll »T a q L . 

(22) 

Irreversible Higgsino annihilation occurs in the HS-phase and governed by the La- 
grangians (|2Tj) and (j2^|) if Higgsino mass is 3> T eW . We calculated all cross sections of 
the t- and s-channel Higgsino annihilation into gauge bosons and massless fermions. The 
technology of these calculations is analogous to calculations in the QCD. The difference is 
we must take into account the annihilation channels, in which initial and final states are 
arbitrary over the two-dimensional color, corresponding to broken SU(2). Full list of all 
annihilation channels is 

XX - BB, XX -> W a W a , 

(23) 

XX -»• B* -> l L l L , e R e R , q L q L , u R u R , d R d R ; XX -> W* -> l L k, qLqL, 

where II, g^, e R , u R , d R is chiral fermions and quarks of three generations. The calculation 
of the total kinetic cross section over all channels (J23|) gives: 

21gj + Qg\gl + 39gl 

M-nn = (24) 

Here g\ = gi(2M x ), gi = #2(2M X ) - gauge coupling constants at the scale y/s = 2M X . 
Their may be calculated by methods of renormalization group analysis. 

From comparison of the theoretical expression for density ()18j) in the assumption, that 
formation of the residual neutralino concentration occurs in the HS-phase, with well known 
experimental constraints (fTTI|) follows that separated neutralino can constitute the Dark 
Matter in the Universe if its mass 

M x = 2900 ± 500 GeV. (25) 

With such mass value the irreversible annihilation starts at the temperature To ~ M x /3 ~ 
1000 GeV, and finishes at the temperature T\ ~ M x /20 ~ 100 GeV . Hence the whole pro- 
cess of the annihilation occurs in the HS-phase that demonstrates the internal consistency 
of the neutralino mass evaluation. 
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Of course, we should treat carefully to the value (|2*5|l since nobody knows what the 
Dark Matter is in reality. Besides this value supposes the execution of M x <C M SUSY . 
However, we should emphasize that the discussing variant of the theory no less motivated 
phenomenologically than well known scenarios with Msusy ~ 1 — 2 TeV. 



6 Problem of separated neutralino searches. 

Above variant with separated Higgsinos can't be tested at LHC because their mass splitting 
is drastically small. So we cannot see any lepton decays of chargino and neutralino. A 
check-up of our model is possible now only in astrophysics. 



6.1 Direct detection. 

6.1.1 Elastic scattering of separated neutralino on nucleons and nuclei. 

There are approximately twenty experimental program for relic WIMPs direct detection 
[2*U] . |21j . [22], [2S1- But we haven't any clear evidence of their existence now. 

All operating detectors use the registration of nuclear recoils. They are generated in the 
processes of the elastic WIMP scattering by nuclei [21]. As usual, the main problem is to 
discriminate signal and background. Now there is a project to use cryogenic apparatus with 
mass near 1 ton in the NUSEL [2E], j2Hj- |2Z|, where the muon background is negligibly 
small. These experiments are planned in the same period that the LHC. 

The most perspective technology is using the liquid xenon as scintillator which react to 
passing nuclear recoils. To study the possibility of neutralino detection in such experiments 
we consider neutralino-nucleon elastic scattering. It is known that spin independent part 
of the total cross section induced by neutralino-quark interaction through scalar quark 
exchange [21j. In our model this part of cross section is strongly damped by large scalar 
quark masses. As a result in the cross section we have only spin dependent part in the 
form: 

9\ml ^(1-4 sin 2 MX 

° xn 64vrm^' °™ Unm^ ' 1 ' 

As it is known the local DM density is 0.3 GeV/sm 3 , the velocity of neutralino at the 
Sun position is near 200 km/ s. Using these data we estimated the separated neutralino flux 
in detector j x ~ 2 • 10 3 sm" 2 s -1 . As example we considered one of the most perspective 
detectors XENON which is in the development stage. The detector threshold is < 10 KeV 
for nuclear recoil energy. The typical recoil energy from interaction with galactic separated 
neutralino is ~ 100 KeV. Detector contains ~ 0.34 m 3 of liquid xenon A = 131.3. 
We calculated the elastic cross section of nonrelalivistic neutralino on nuclei of xenon ~ 
10~ 35 sm 2 . The estimation of number of events gives the value that close to the planned 
background for considered detector 4.5 ■ 10~ 5 s _1 or ~ 1 event per hour that can't provide 
the reliable signal for existing of Dark Matter in the form of separated neutralino. The 
question about direct detection possibility of separated neutralino remains open. 

6.1.2 Recharging of separated neutralino on nucleons. 

New and exotic processes could be possible when relic neutralino penetrate through the 
matter. Corresponding channels are: 
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Cross sections have the form: 



g\m 2 v jE N — Am N — AM X 



a 



x '' W~n,( v V /'. \ 

,4^,2 



#2 m p E N + Am N -AM : 



x 



xn ldnm^ V 

Here Am^v = — m p — IMeV - proton-neutron mass difference; AM X = M^± — 
M x o ~ 360 MeV" - chargino-neutralino mass difference that can be calculated from ()13j) ; 
En = \rri]s[vl el . An average kinetic energy of neutralino in the locality of the Sun is about 
~ 1 MeV, but recharging reaction is possible only if there are high energy neutralinos with 
Eki n > 1 TeV in cosmic rays. 



6.2 Indirect detection. 

6.2.1 Diffuse gamma ray spectrum. 

Except of direct relic neutralino detection experiments there are many projects use satellite 
experiments for the neutralino annihilation products observation [20], [22], [21] • It can be 
possible to establish some qualitative and quantitative WIMP's characteristics from these 
data, mainly from the data on the neutralino annihilation spectrum from the Galactic 
halo. We consider this process and calculate characteristic energies in the diffuse gamma 
spectrum. Dirac neutralinos in the halo annihilate into W- and Z-bosons 

XX - W+W-, XX - ZZ, (27) 

and into fermions 

XX -> qq, ll (28) 
Corresponding total kinetic cross section given in the form: 

£ 2 4 (21 - 40cos 2 fl w + 34cos 4 M 
[(TV)ann 256 ttMW^ 1 } 

It is known that total multiplicity of secondary hadrons nearly twice larger than charged 
hadrons multiplicity |2Hj- [SOI- Average multiplicity of charged secondary hadrons 
(n-ch) (V~s) was studied in e + e~ , pp and pp reactions, e ± p interactions . It was established 
that charged multiplicity is an universal function of energy 

n>ch{Vs) = A + B ln + C In 2 y/s, n ch = (n ch ) (y/s / q ) - n 

with fixed parameters 

A = 3.11 ± 0.08, B = -0.49 ± 0.09, C = 0.98 ± 0.02. 

To choice a specific channel it is need to fix parameters qo, n . For neutralino annihilation 

Voted = X > n o(xx) = °- 

We have used experimental data on multiplicity in Z — > hh decay For estimation 
we suppose that charged hadrons part x = (n ch ) / (n h ) doesn't depend on energy and have 
the value x ~ 0.49. 
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To plan an experiment it is important to know characteristic photon energies that are 
generated in the following decays: 7r° — > 27, if — > 37, rf — > 3tt° — > 67. In the annihilation 
channel (}28|) the total hadron multiplicity is described by logarithmic function well: 



(n f h f ) = x~\A + Bln2M x + Cln 2 2M x ) ~ 149, M x ~ 3 TeV. (30) 

Average energy of the neutral pion in neutralino annihilation products is E n o ~ jE^o ~ 
2M x /{n(/). So we have characteristic maximal photon energy in the n° — > 27 decay: 

E 7{ ^ 27) ~ E n o/2 = 20 GW. (31) 

Analogously from decays 77 — > 37 and 77 — > 37r° — > 67 energies of photons are: 

J B 7(r) o_3 7) ~ ^0 /3 = 13.5 GeV , £ 7(7? o-6 7 ) ~ E v o /6 = 6.7 CeV. (32) 

Then we considered second annihilation channel (|27|). Total hadron multiplicity of W- 
and Z-bosons decays is 

(nY Z ) ^ 42.9. (33) 

Average energy of the neutral pion is E w o ~ i^o ~ M x / (nf[ z ) . So we have characteristic 
maximal photon energies: 

£ ! 7 ( 7r o_ > 2 7 ) — E w o/2 ~ 35 GW; 

(34) 

£ 7( „o^3 7) ~ ^0 /3 = 23.3 GeV , £ 7(7? o^6 7 ) ~ E v o /6 = 12 GeV. 

Multiplicity distribution in a wide energy region is described by the negative binomial 
distribution (NBD) with a good precision, which depends on energy very weak (as loga- 
rithm) |H]: 

, _ fc(fc + l)...(fc + n- 1) (h/k) n 
(n;n ' } ~ n\ ' [l + {n/k)Y+ k] (35) 

k' 1 (a/s) = a + 6 In 

where n = n c h, n = n c /j. Coefficients of the function A;~ 1 (y / i) is different for various 
channels: 

a e+e - = -0.064 ± 0.003, b e+e ~ = 0.023 ± 0.002; (36) 
a P p/pp = -0.104 ± 0.004, b pp/Pp = 0.058 ± 0.001. (37) 
Using NBD we found approximate photon distribution over the energy. The number of 
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photons with the energy E 1 in one xx-annihilation act may be defined by 
dAL 2M V 



dE~ 



{ [Br(7T°/h) + Br(rf/h) ■ Sr(r/° -> 2 7 )] x 



Br(h)(nH)P ( ^ (njf),*,,) + l -Br{WZ)(n^)P f^- 



( n ch Z )i k 



wz 



+Br(rf/h) 



Br{rf -> 3tt°) + -5r(^ tt+tt-tt ) 
3 



X 



Br{h)(n{l)P ( ^ <*#>,*//) + ^r(^)(nDP (-^; « z ),^ 



x 



+ -Br(7]°/h)Br{r]° -> tt+tt^x 
3 



M. 



6£. 



+ 



+ 



(38) 

Here Br(WZ) ~ 0.2 is the total branching for neutralino annihilation into W 7 - Z-bosons; 
charge hadron multiplicities and {n^ h z ) defined in f)3Up and values /cjj 1 = 

k~ 1 (2M x ) = 0.4 and fc^ = k~ l (M x ) = 0.12 defined with coefficients flH7|) and ® corre- 
spondingly. The spectrum is shown graphically at Fig. El 



dN/dE, 
[1/GeV] 




30 40 
E, [GeV] 

Figure 2: Diffuse gamma ray annihilation spectrum. 



Calculated gamma-ray flux of annihilation radiation from clumped galactic halo is closes 
to sensitivity threshold of the satellite detector GLAST. The question about principal 
possibility for registration of the neutralino annihilation spectrum at GLAST remains open. 
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6.2.2 High energy gamma ray bursts as a result of collapse and total annihi- 
lation of neutralino clumps. 

Consideration neutralino as a carrier of DM mass in the Universe generates questions about 
their evolution and about astrophysical phenomenons which evidence both the existence of 
neutralino and the reality of the evolution process with neutralino participation. 

According to high resolution cold dark matter simulations in framework of hierarchical 
principle of structure formation large virilized halos are formed through the constant merg- 
ing of smaller halos formed at earlier times. Neutralinos at the gravitational isolation stage 
z = 8 -T- 10 formed dense massive clumps of different scales j3H|, [HI]- It likes naturally that 
DM clumps are in the state of dynamical equilibrium as a result of compensation of own 
gravitational forces by forces of tidal interactions with each other. In consequence of inner 
instability the clump can turn into irreversible collapse state. The growth of density pro- 
vides the growth of annihilation rate. So neutralino clouds become more and more intense 
sources of relativistic particles. At final stage of this process there is a gamma ray burst 
with maximal characteristic photon energies defined above (JHH), (JH2I), (J33|) . Specific energy 
and time distribution of the gamma radiation are connected with neutralino parameters 
and give an important information about structure of the MSSM and about its correspon- 
dence with cosmology. Therefore our main statement consist in the nature of rare powerful 
cosmological gamma ray bursts is a result of the total annihilation of neutralino clumps. 

We analyzed the qualitative and quantitative descriptions of this process in the spheri- 
cally symmetric collapse model with two assumptions. First, annihilation products leave a 
clump during a time substantially smaller than the time of its macroscopic evolution. And 
second, an annihilating clump is spatially homogeneous and isotropic since more dense 
regions annihilate faster than the rest and heterogeneities are vanishing. 

The model, defined by equations of mass and momentum balance and the equation of 
nonrelativistic gravitational theory, is: 

- J Px dV = -J Px v ld S t - J ^ dV , 

^ j Px v idV = - j {p x ViV k + P ik )dS k - j p x Vi<pdV - j v dis p x VidV - j ^ nnPx v id V 

A0 = 47rG'p x , 



(39) 

where Pi k is the pressure tensor for neutralino gas. It is absent in consequence of second 
assumption. Last terms in balance equations describe the momentary going away of an- 
nihilation products from the collapsing clump. Taken into account that annihilation and 
scattering are cross-channels. In the local form we have the simple system: 

dM 3M m „ M 2 fdR\ 2 2GM / „ 9 R 



P 2 jr , ( 4 0) 



dt 2nM x # 3 ' \dt J R V #o 

which integrates elementary. Here R is a clump radius; M = 4irp x R 3 /3 is a clump mass; 
Ro = const is an initial radius; (3 2 < 1 is the constant which fixed by initial velocity of the 
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surface clump compression. The results are in the parametric form 



M(C) = M W 



_ {(Tv)ann ( 3p X (*) \ ^ ,2 

M x { 2nG J V 3 " 



-2 



R(0=Ro(( 2 + f3 2 r 1 



, Mann / 3p X (*A ^ A ((2 + ^ + 2/? 2 (( 2 + F ~ 1) 



3M X V J \ ^ ^ ' C 2 + P 2 

(41) 

where M(#) ~ M = 4 • 10 33 <7, p x (*) ~ 10 4 g/cm? are initial clump mass and density. As 
you see M(t) -> 0, i2(t) -> for i -> oo. 
Annihilation rate is 



»/ 2M awn (C 2 + /3 2 ) 3 

= ^ann - ^(*)P X (*) 17 7 77^ 74" ^ 42 ^ 



Generally behavior of this function is that a slow growth at small times changes into quick 
evolution around narrow peak after which there is an exponential drop to zero. This 
behavior explains by changing of the parameter hierarchy during the collapse. We use the 
convenient time scale 

• _ M(„) / (t-t ) 2 \ _ Mann . 

Eann - — \ l -^T)' T °-^GM~ X (43) 
Here to is the time of maximum. As we have 

= (Harm ( 3p X (*) \ ^ ~ , 

* M x V 2vrG / 

than the annihilation rate doesn't depend on initial conditions and expresses through 
fundamental constants only. Dependence of function E ann jMu\ via dimensionless time 
t] — (t — t Q )/r shown at the left Fig. El Special features of this function is skewness. 
About 80 % of clump mass annihilates during the time 



2(cxu) ann 0|(21 - 10 cos 2 6 W + 34 cos 4 6 



t ann = 6r = = g ^ io — ^ 30 s , M x ^3Tey (44) 



7rGM x 128 ttWGcos 4 ^ 



at that effectively it proceeds between = Atq (before maximum) and = 2tq (after 
maximum) . 

The mass evolution M x is shown at the right Fig. |3J Decay time of collapsing neutralino 
clump expresses through MSSM parameters and gravitational constant but it is absolutely 
independent from initial density and initial velocity of clump compressing. 
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Figure 3: Annihilation rate Ej^/Mm and relative mass M/Mr^ of neutralino clump as a function 
of dimensionless time rj. 

7 Conclusion. 

Here we discuss the variant of supersymmetric extension of the SM, where Higgsino mass 
scale is lower than the scale of soft SUSY breaking: // <C M S usy ~ M ~ Mi/ 2 - We call 
the model as Model with Separated Higgsinos. 

In the model, as in the commonly used MSSM, the precise convergence of invariant 
charges takes place at the scale that does not contradict to proton mean life time con- 
straint. When Msusy > 10 7 GeV radiation corrections dominate in the neutralino-chargino 
mass splitting. From the calculation of splitting and the renormgroup evolution we get 
M x ± — M x o 2 ~ 360 MeV at M x ~ fi « 3 TeV. This value have got from the analysis of the 
neutralino relic in cosmological plasma. We supposed that neutralino is a DM carrier. In 
the model formation of the neutralino relic occurs in high symmetric phase of the plasma. 
Charged Higgsino mean life time is r ± ~ 0.4 x 10~ 9 s. Problems of relic separated neu- 
tralino search in underground (NUSEL) and satellite (GLAST) experiments are discussed. 
It have been shown that spin dependent component dominates in neutralino-nucleon cross 
section. Possibility of the separated neutralino registration at perspective detector XENON 
is analyzed. In this case the signal is closed to the planned background. Process of high 
energy neutralino recharging on nuclei have been considered for the case when neutralinos 
with large kinetic energies (more than 1 TeV) are in cosmic rays. The spectrum of diffuse 
annihilation radiation have been calculated and analyzed for the galactic DM. Annihilation 
cross sections for galactic neutralino have been calculated. Contribution to the gamma ray 
spectrum from t- and s-channel neutralino annihilation into Z°Z° and W + W~ pairs has 
been calculated using experimental data about hadron multiplicities at mass surfaces of 
Z°, H^-bosons; quark- ant iquark s-channel contribution has been estimated using the phe- 
nomeno logical model of hadron multiplicities at ^Js = 2M X . In the Galaxy a collapse and 
total annihilation of neutralino clumps should take place if neutralinos are main carriers 
of its mass. Time and energy distribution of the gamma ray radiation are calculated. It is 
shown that nearly 80 percents of the clump mass annihilate in 30 seconds. 
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